CALT-68-2058 

DOE RESEARCH AND 

DEVELOPMENT REPORT 



Mean Field Theory for Fermion-based U(2) Anyons 



Patrick McGraw^ 



California Institute of Technology, Pasadena, CA 91125 



Abstract 

The energy density is computed for a U(2) Chern-Simons theory coupled to a non- 
relativistic fermion field (a theory of "non-Abelian anyons") under the assumptions 
of uniform charge and matter density. When the matter field is a spinless fermion, we 
find that this energy is independent of the two Chern-Simons coupling constants and 
is minimized when the non-Abelian charge density is zero. This suggests that there 
is no spontaneous breaking of the SU(2) subgroup of the theory's symmetry, at least 
in this mean field approximation. For spin-1/2 fermions, we find evidence of ground 
states with a small non-Abelian charge density, which vanishes as the theory of free 
fermions is approached. 
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1 Introduction 



In this paper, we consider a 2 + 1 dimensional model with a non-relativistic matter 
field \l> minimally coupled to an SU(2) x U(l) — U{2) "statistical" or Chern-Simons 
gauge field. The matter will be taken to carry a unit of U(l) charge and a fundamental 
representation of SU (2) "isospin." The Lagrangian is given by: 

(1) 

where and A a ^ are the U(l) and SU(2) gauge fields, with Chern-Simons coupling 
constants k and k respectively, and the covariant derivative is 

D^d^ + iA^ + iA"^. (2) 

(It will sometimes also be convenient to define = A^ + A a ^a a /2.) We will consider 
cases of either bosonic or fermionic matter fields, so that \l/ may obey either canonical 
commutation or anticommutation relations. The Hamiltonian corresponding to (1) is 
given by 

H = J d 2 ^(D^(D^) (3) 

while the gauge fields are subject to constraints which relate them to the matter 
fields: 

B = -F u = eijdiAj = --p, 

\ a i 1 abc \ b ac\ ■ 



B a = = eijidiA" + -e abc A\A^ = -^p\ (4) 



and 



Foi = --eijJj, 

K 



Foi — ~^ e ijJj^ (5) 
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where the densities and currents are defined by: 




Since there is no Maxwell term, the gauge fields are non-dynamical and are completely 
determined (up to gauge transformation) by (4) and (5). 

As usual, the theory is gauge invariant only when k is an integer multiple of 1/Att. 
The Gauss constraints ensure that any electrically charged particle also carries mag- 
netic flux, so the particles experience mutual Aharonov-Bohm interactions resulting 
in exotic statistics. In this case, the braiding of two particles leads not only to an 
overall phase change of the wave function but also a non-Abelian rotation acting on 
the isospin indices of the two particles. A two-particle total wave function ip(xi,X2) 
is acted on by the braiding operator 



where can be thought of as a matrix-valued SU (2) "phase." The eigenstates of this 
two-particle braiding operator are states of definite total isospin, i.e., pure iso-singlet 
or triplet states. 

This model of "non-Abelian anyons" was studied in reference [1] using a mean field 
approximation. In this mean field technique one first quantizes the matter field in the 
presence of a classical background magnetic field. One searches for a self-consistent 
ground state having a uniform expectation value of the matter and isospin densities 
consistent with the Gauss constraints. Then, using the densities and currents as fun- 
damental variables, one can study the effect of fluctuations about this approximate 
ground state, using the Bogoliubov approximation in which the commutators of fluc- 
tuation operators are replaced by their expectation values in the zeroth-order ground 




(7) 



state. 
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In reference [1], it was shown that with bosonic matter the mean field ground 
state energy behaves differently depending on the values of the two Chern-Simons 
couplings. In particular, there was a phase in which the energy was minimized in 
the mean field approximation by the generation of a non-zero isospin density, and 
thus the SU (2) symmetry was spontaneously broken. This paper will investigate the 
consequences of coupling fermions instead of bosons to the U (2) Chern-Simons field. 
We will show that for spinless (or polarized) fermions, the resulting mean field theory 
differs from the bosonic case in two ways: (1) The form of the mean field energy 
density as a function of matter and charge densities is independent of the Chern- 
Simons coupling constants and (2) For a given particle density, the lowest energy 
always occurs when the SU(2) charge density is zero. Thus the mean field picture 
does not show any breaking of the £77(2) symmetry. In the case of spin-1/2 fermions, 
however, there is a hint of spontaneous symmetry breaking which vanishes smoothly 
as the theory approaches that of pure fermions. 

2 Mean field approximation for fermion-based U(2) anyons 

For a theory where the matter field \& is fermionic, we now search, much as in 
reference [1], for a self-consistent ground state \ fl) with uniform matter and isospin 
densities: 

(p) = PO, (Ji) = 0, 

(P a ) = Pi (J?) = 0. (8) 

According to the Gauss constraints, these nonzero densities will imply uniform mag- 
netic fields. We will consider matter in the background of these fields, and then 
demand that the number of particles per unit area be consistent with the Gauss con- 
straints. We can take the isospin density to be along the 03 direction, (p a ) = 5fpo- 
Then we have U(l) and £77(2) magnetic fields given by 

(B) =B = — -po, (B a ) = 5 a s B = -l%po. (9) 

K K 
[3] 

In the symmetric (or isotropic) gauge A% = (B/2)eijX J , the gauge field can be 



4 



written as 



A = e 4jT (- + - ( x 3 ). 



(10) 



The single-particle orbitals are split into two sets according to the eigenvalue of 
03. Particles in states of up or down isospin feel different effective magnetic fields, 
B+ = (Bq + Bq/2) and B- = (Bo — Bo/2), respectively. The single-particle energy 
levels thus fall into two sets of Landau levels. The kth energy levels of the isospin up 
and down systems have energies given by 

* = 5). in) 

and have degeneracy per unit area N±/A = \B±\/2n. 

Uniformity of the density requires that all orbitals of any given Landau level be 
filled with the same number of particles.' 41 In the case of spinless (or spin-polarized) 
fermions, this means either or 1 per orbital. For spin-1/2 fermions, the possible 
occupancies are 0,1, and 2. Let us first consider spinless fermions. The total energy 
per unit area for a state with the lowest n levels filled is given by: 

A ~ 2tt ^ Cfc ~ A-nm ' ^ ' 

(The factor B /2n in front of the sum represents the degeneracy.) In our system, if 
the lowest n + and n_ of the isospin up and down levels, respectively, are uniformly 
filled with one particle per orbital, then the combined energy is given by 

E - « 2 + (Bo + k? + i!l (Bo _ 4,2. (13) 



A 4vrm v 2 ' Anm K 2 

The density of isospin-up particles is given by the number of filled Landau levels 
times the degeneracy, p+ = n + N + /A, and similarly for the down particles it is p_ = 
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n-N-/A. The matter density p is the sum of the up and down densities, while the 
isospin density for isospin-1/2 particles is given by half the difference' 11 : 

n+N+ n_iV_ \B + B /2\ \B - B /2\ 

P0 = p+ + p - = -4- + -4- = n+ 2, + n - Y ' 

- l ( , 1/ \By + B /2\ \Bo-Bo/2\ 

Po = ~ 2 {P + ~ P-) = ~ 2 {n + n_ ). (14) 

These equations may be rewritten in the form 

n+\B + B /2\ po _ n_\Bo-B Q /2\ p _ . . 
Y = T + P °' 2~n = T - P °- (15) 

The expression (13) for the energy then becomes: 

For comparison, the result in the bosonic case of [1] was found to be: 



E 1 

A 2m 



i Po . po w po , ~ x . i PO PO ,/PO ~ x, 
|- + ^l(y+Po) + |--^l( y -Po) 



(17) 



The latter expression is different because in the bosonic case, only the lowest of 
each of the two sets of Landau levels is occupied in the ground state, whereas in the 
fermionic system, the exclusion principle requires that higher levels be occupied. The 
two expressions (16) and (17) are plotted in figure 1 for representative values of the 
coupling constants. 

We note the following features of the fermionic result which differ from those 
of the bosonic result: (1) The expression is completely independent of the Chern- 
Simons coupling constants k and k, and was derived without any reference to the 
Gauss constraints —kBq = po and —kB$ = po. (2) For a given matter density po, 
the energy is always minimized by po = 0. Thus it appears that the assumptions 
(14) for the ground state are not self-consistent unless po = 0, and there is, in this 
approximation, no spontaneous breaking of the SU (2) symmetry. 
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Figure 1: Upper plot (A) (adapted from reference [ 1]) : energy of bosonic 
mean-field theory plotted against the ratio po/po f° r three values of the coupling 
constants. Curve I is for k/k — 0.5, and is typical of the region kk > 0. Curve II, 
with k/k = —0.1, is typical for kk < and < .25. Curve III, with k/k = —0.7, 

is representative of kk < 0, > 0.25. On the x-axis is the ratio po/po, on the 

y-axis is energy in units of l/2m/t. In case I, The minimum occurs at po = 0. 
In case III there is a self-consistent mean-field ground state with maximal isospin 
density, po = po/2 In case II, the mean-field approximation is not self-consistent: The 
minimum energy occurs in a limit where one of the effective magnetic fields B+ or 
B- goes to zero and the Landau level picture breaks down. The lower graph (B) 
shows the much simpler form of the energy expression for spinless fermions: it is 
independent of k and k, and is always a minimum when po = 0. The energy scale for 
the lower plot is ii/2m. 
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The first of the above observations is not entirely surprising in view of the results 
for the Abelian model (which should correspond to the limit <C 1/I K D- When 

po = 0, (16) reduces to the result for fermion-based Abelian anyons. This result is 
similarly independent of the Chern-Simons coupling, and is equal to the energy of a 
degenerate fermi gas in 2 + 1 dimensions. Corrections to this fermi energy are found 
only when one includes the effects of fluctuations to quadratic order. The second 
property, the absence of a spontaneous non-Abelian charge density, is less obvious. 
In the boson-based case, it was found that the phase with po = ±po/2 was stable 
in the regime l/\k\ < 4/|/c| with k and k having opposite signs. Fermions can in 
principle be generated from bosons by setting the Abelian Chern-Simons coupling 
to k — —l/2n, resulting in a statistical angle of tt. If we were to continue naively 
from the behavior near the bosonic point (k — > oo) to k — —1/2tt, we might expect 
symmetry breaking at the fermionic end when \k\ > g^. It might be argued, in view 
of the spin-statistics connection, that it is not natural to expect a theory of spinless 
bosons to be connected continuously to one of spinless fermions. The properties of 
the pure Abelian theory obtained with the prescription of refs. [1] and [2] interpolate 
smoothly between the bosonic theory and a theory of unpolarized spin-1/2 fermions.' 51 
Therefore, in the next section, we will consider a spin-1/2 matter field. 

Qualitatively, the comparatively greater susceptibility of the bosonic model to the 
development of a spontaneous asymmetry between the isospin-up and down Landau 
levels may be explained in terms of the particles' exclusion properties. An isospin 
asymmetry generally raises the energies of one set of Landau levels while lowering 
those of the other. Since any number of bosons may occupy the lowest Landau 
level, it can become energetically favorable to lower the energy of, say, the isospin-up 
Landau level and place all of the particles into this lowered level, while leaving the 
raised level unoccupied. The fermionic system, on the other hand, can be regarded 
to a first approximation as two fermi fluids, one filling the isospin-up levels and one 
filling the down levels. The development of an isospin asymmetry requires that more 
particles be added to one of the two fluids. Even if energy of the kth Landau level is 
lowered by the asymmetry, the levels must be filled up to a higher value of k, which 



8 



offsets the energetic advantage. 



3 Spin- 1/2 fermions 

We now ask whether the result of the previous section is changed if we consider 
spin- 1/2 fermions instead of spinless ones. For spin- 1/2 fermions, there are two states 
per orbital, and so it is possible to fill a Landau level with either 1 or 2 particles 
per orbital. (We are supposing that the statistical gauge field does not couple to the 
spin, so that this double occupation is the only effect.) Let n — 2p + <r, a — 0, 1, and 
consider a state in which the lowest n/2 of a set of Landau levels are filled. If n is 
odd (a = l), we mean by this that the lowest p levels are doubly filled, and the p + 1 
level is filled with one particle per orbital. The total energy per unit area of such a 
configuration is' 11 

E B , v — -\ B , , 1. B , 1.. B 2 , o . .. 

— = —\2y—(k + -)+a—(p + -)] = (n 2 + a . 18 

A 2tt ^ rn 2 J wT 2 n 8nm K J K ' 

For our system with two sets of Landau levels, the total energy becomes 

E = (nl + , + ) Bo 2 (nl+„-) B» 

A 8vrm V U 2 ' 8nm K 2 J K ' 

As before, the + and — subscripts refer to the isospin states. Whereas the corre- 
sponding equation (13) for the spinless case only involved the products u±\Bq ± -^1 
and thus could be expressed in terms of the densities without reference to the Gauss 
constraints, that is not the case here. Using the Gauss constraints, we write 

Bo± T = ~ T 2k> 



and thus 



E _ n \(Pl , ^ , a + ,P0 , P0 x2 , °- ,P0 P0 n2 



Since the second and third terms are non-negative, we see that, at the mean-field 
level, the lowest energy state for a given po is still one with a + = c_ = and po = 0. 
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However, if a + or o"_ is restricted to be 1, then a minimum of the energy does in fact 
appear at 

PO = ± — 2 - 2 , J PO ~ 1fl 2 ~ ' 21 

K 327T Z K Z + 1 Ioh^kk, 

(See figure 2.) This expectation value of p vanishes in the limit of free fermions 
— > 0, 1/k— > 0). The corresponding energy, in the limit where both and 1/k 
are small, is given by 

!^7 + 8^ + °<> < 22 > 

It is worth noting that in the Abelian case, states with the top level half-filled are 
also energetically unfavorable at the mean-field level, but that fluctuations introduce 
a- dependent corrections of the same approximate size (0(1/ k 2 )) and opposite sign. 
Thus it is conceivable that in our model, the a± — 1 ground states might be stabilized 
by quadratic corrections. Also, in the Abelian model, the half-filled ground states are 
the only consistent ones at odd values of the coupling constant. 



4 Consistency of mean-field theories 



The assumption of a uniform matter density in a uniform magnetic field requires 
that all orbitals within a Landau level be filled equally. This means that the only 
degree of freedom for the ground-state distribution of particles in a set of Landau levels 
is the filling factor, which must be an integer. In the case of a bosonic ground state, 
this integer represents the occupation number of all orbitals in the lowest Landau 
level, while in a fermionic ground state it represents the number of levels which are 
filled. Combined with the Gauss constraints, this condition will in general pick out 
a discrete set of values of the Chern-Simons coupling constants at which the mean- 
field ground state is well-defined. In the case of Abelian anyons (either boson- or 
fermion-based) there is only one set of Landau levels which may be filled with an 
integer factor. The matter density is given by the filling factor times the Landau 
level degeneracy per unit area. The degeneracy is related to the magnetic field, which 
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Figure 2 A: Energy of fermion- based anyons for a + = <r_ = (top) and for 
<r+ = 1,(T_ = (bottom). 
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Figure 2 B: The two curves of figure 2 A, shown on the same axes. 
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is in turn related to the density by the Gauss constraint: 

n \ B o\ n po 
P0 = — = ^- (23) 

po m ay be divided out from both sides, leading to the familiar series of allowed values 
of the coupling constant: k = n/2n, corresponding to the series of statistical angles 
2ir/n and 1 — 2n/n based on bosons and fermions, respectively. 

Similar constraints occur in the non-Abelian model we are considering, but the 
consistency conditions are more complicated because two sets of Landau levels are 
involved. This issue was not addressed in reference [1], so we consider here both the 
bosonic and fermionic systems. The basic equations are those of (14), which relate 
the isospin and matter densities to the filling factors and magnetic fields. When the 
Gauss constraints (4) are applied to substitute for the magnetic fields in (14), these 
equations become: 



n + i Po , Po i , n-.po po 



2n k 2k 2h k 2k 



Ib +,H0 l J I '<<- ,/AJ HO I / .n 
P0 = tH— + 7d-7H— -7^|- ( 24 ) 



These basic equations have the same form for fermi and bose-based systems; only the 
interpretation of the filling factors n + and n_ is different. 

Now consider a mean-field state with 

V 

Po = 2^0, 

where — 1 < rj < 1. In such a state, equations (24) become: 

1 n+ ^ + ^ + n ~ ^ ^ I 
2n k Ak 2ix k 4k 

V = ^ ± \ 1 - + ^\--\--^ L \. (25) 

Because of the quantization of k we may write k = m/Air, where m is an integer. 
With this substitution, the equations (25) may be combined into a relation between 
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n + and n_, and another relation involving only n + : 

d-^i- + ^ = (i+^i--^i, (26) 

2n re m 2n k m 

^ = |- + ^|. (27) 
n + rem 

In the case rj = 0, these reduce to k = y = y, reproducing the familiar set of 
Abelian mean-field theories. The bosonic ground state of maximal isospin alignment 
described in ref. [1] corresponds to rj — 1. In this case, we have n_ = and 

I = ^-i, (28) 
re n + m 

and thus find that re takes values which are rational, but not necessarily integer, 
multiples of re approaches integer values only in the limit m ^> n+, or re ^> re. 

We now consider the mean-field energy minimum (21) of the spin-1/2 fermion- 
based theory, which corresponds to 

-4re / 1 \ -m ( 1 \ 
V ~ ~V V32vr 2 K 2 + I )--^ [ 2 m 2 + 1 J ' { } 

The conditions on the coupling constants for a consistent mean-field theory at this 
value of i] turn out to be more complicated. Noting the useful expressions 

' 7rre(2m 2 + l) V ; 



and 



1 77 . 1 / 2m 2 + 1 =F 1 



re ± 4k' Irel 1 2m 2 + 1 ' ' ( ' >j) 



we find that equation (26) (which relates n_ and n+) takes the form 

n + [7rre(2m 2 + 1) + m]m 2 = n_[vrre(2m 2 + 1) - m](m 2 + 1). (32) 

Note that, as usual, 7rre must be rational; we may write 7rre = p/q, where p and q are 
relatively prime integers. By assumption, n_ is even while n + is odd. When equation 
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(32) is multiplied by q, the RHS is even, so that n + [p(2m 2 + l)+qm]m 2 must likewise 
be even. This can only be satisfied if either m is even, or p,q, and m are all odd. 

If k > 0, the other consistency equation (27) becomes, after substitution of the 
expressions (30) and (31): 

7TK(2m 2 + 1) — m — 2m 2 n_, 



or 

2m 2 n+ + m 

7TK = o • (33) 

2m 2 + 1 K ' 

Substitution of the above expression for ttk in (32) yields: 

m(mn + + l)=n_(m 2 + l). (34) 

Letting n + = n_ + D, with D an odd integer, we find 

D{m 2 + 1) + m = n+. (35) 

By substituting this into the expression (33) for ttk, we finally obtain a relation 
between the two coupling constants: 

2m 2 ((Dm 2 + 1) + m) 

nK = n — 9 1 

2m 2 + 1 

Dm 2 = l6ir 2 Dk 2 , (36) 



the latter expression being valid in the limit of small inverse coupling constants ^, A 
(near the point of free fermions). 
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5 Discussion and Conclusions 



The relation (36) describing the coupling constants at which ground states with 
broken SU(2) occur is a rather peculiar one. It is nonlinear, and does not appear 
to connect continuously to the theory's behavior near the bosonic point. Translated 
into the variables 9 = 1/ttk and 9 = I /irk, which represent the sizes of the Abelian 
and non- Abelian statistical phases, the relation becomes (in the 9, 9 — > limit): 

9 = — -. (37) 
16D K ' 

In figure 3, the behavior of the mean field theory is plotted in the 9 — 9 plane 
between 9 = and 9 = it. The points (0, 0) and (it, 0) mark free bosons and free 
fermions, respectively, and it is near these two points that one expects the mean-field 
approximation to be useful. The odd-n fermion-based ground states appear between 
the line 9 = n and a parabola, and are clustered near the 9 = n axis. There is no 
obvious way to continue this behavior to that at the bosonic point. 

One might draw one of two conclusions: Either there are phase transitions on 
the 9 — 9 plot other than the ones shown (i.e., between 9 = and 9 = ir), or mean- 
field theory alone is not sufficient to understand the symmetry-breaking behavior of 
this theory near the free-fermion point. The fact that the expectation value (21) of 
the isospin density is second order in inverse coupling constants lends credence to 
the suspicion that the second conclusion is true: mean-field results typically receive 
corrections at quadratic order in inverse coupling constants due to fluctuations. The 
results in this paper, however, suggest that: (1) if there is a spontaneous isospin 
density near the fermionic point, it vanishes smoothly as that point is approached, 
and (2) it is worthwhile to study the question using other methods. It has been shown, 
for example, that the pure SU(2) theory (to which our model reduces in the limit 
k/k — > 0) is susceptible to the formation of Cooper pairs in an isosinglet state.' 61 The 
condensation of isosinglet Cooper pairs would naturally form a ground state with zero 
isospin density. But it is possible that when the Abelian coupling is also included, 
isotriplet Cooper pairs might form in some regions of parameter space. 
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Figure 3: The 9 — 9 plane. The theory's behavior near the points (0, 0) (pure 
bosons) and (0, 7r) (pure fermions) delineates several regions. In region I, there is no 
spontaneous breaking of SU(2): the energy is minimized by p = 0. In region III. 
there is a consistent mean-field ground state with a maximal expectation value of 
(p) = Po/2. In region II, mean-field theory fails: The ground state energy appears to 
have a minimum at some (p) < po/2, but no self-consistent state can be constructed 
at that minimum. These are the three regions of parameter space corresponding to 
the curves I, II, and III in figure 1. Near the fermion point (0, 7r) (regions IV and V), 
there are mean- field ground states with (p) = 0. In region IV, there also exist states 
differing from the others only at quadratic order, and having a small expectation 
value of p. These are the odd filling factor states we have been discussing. It is not 
clear how these regions might connect to one another in areas far from the fermionic 
and bosonic points. 
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